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Aims. We present mid-infrared observations and photometry of the transitional disks around the young stellar objects DHTau, 
DM Tau, and GM Aur, obtained with VISIRA'^LT in N band. Our aim is to resolve the inner region and the large-scale structures 
of these transitional disks, carrying potential signatures of intermediate or later stages of disk evolution and ongoing planet formation. 
Methods. We use the simultaneously observed standard-stars as PSF reference to constrain the radial flux profiles of our target objects. 
Subtracting the obtained standard-star profile from the corresponding science object profile yields the flux residuals produced by the 
star-disk system. A detection threshold takes into account the background standard deviation and also the seeing variations during the 
observations to evaluate the significance of these flux residuals. On the basis of a simple model for the dust re-emission, we derive 
constraints on the inner radius of the dust disk. 

Results. We spatially resolve the transitional disk around GMAur and determine an inner-disk hole radius of 20.5^Q g AU. The cir- 
cumstellar disks around DH Tau and DM Tau are not spatially resolved but we are able to constrain the inner-disk hole radius to 
< 15.5^20 and < 15. S;!;"^ AU, respectively. The performed photometry yields fluxes of 178 + 31 mJy for DHTau, 56 + 6mJy for 
DM Tau; and 229 ± 14 mJy for GM Aur 
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1. Introduction 

Classical T Tauri stars are low-mass, pre-main sequence stars 
surrounded by circumstellar disks. These disks are a natural out- 
come of the star formation process and are known to dissipate 
over time by means of stellar winds or photoevaporation by 
either the radiation of the central star or some external radia- 



tion (e.g., Hollenbach et al. 2000[parke et al.|2001, Alexander 



& Armitage 2007 1, accretion onto the central star (Hartmann 
et al.lllQgS*), grain growth and fragmentation (e.g., [van Boekel 



et al.||2004, DuUemond & Dominik 2005), and the formation 



of pl anets (e.g., [Lin & Papaloizou 1986; Marsh & Mahoney 



1992 



et al. 



Pollack etal.[T996, Bryden etal..l999, Boss, 2 002 ; Calve| 



neyl 
I vet] 



bris disks found around main-sequence stars (e.g., |Calvet et al. 
'20051 |DWessio|2009| ). The spectral energy distribution (SED) 
of transitional disks shows a deficit of emission in the near- 
infrared, steep slopes in the mid-infrared and an excess at far- 
infrared wavelengths that can be explained by inner-disk gap(s) 
or hole(s), while a prominent outer disk is still present (e.g., 
[D'Alessio et al.|2005| ). To date, the structure of the transitional 
disks has been mainly investigated by model fitting of their 
SED, which usually does not provide an unequivocal solution 
(D'Alessio et al. 2005 ; Calvet et al. 2005 ; Ireland & Kraus 2008) 
jD'Alessio 2009) . Imaging of protoplanetary disks provides an 
independent check of SED modeling and a far less ambiguous 
method for studying the disk structure (e.g ., [Chen & Jura|2003 



important implications for the possibility of planet formation. 
The key to constraining the physical mechanisms responsible for 
disk evolution is the detection and detailed analysis of disks that 
are in the process of dispersing. The evolution of circumstellar 
disks is mostly studied in the infrared by analyzing stars that 
exhibit emission in excess of that expected from the stellar pho- 
tosphere. In contrast to the case for high-resolution imaging in 
the optical/near-infrared, the contribution from the stellar pho- 
tosphere is small compared to the disk emission at mid-infrared 
wavelengths. 

Transitional disks are considered to represent an evolved 



2002]). The way m which the disk material dissipates has i Pantin et al.|2005, Eisner et al.. 2009, Hughes et al.|2009t|Sauter 



stage of the accretion disks s urrounding a Class II source (Lada 
|& WilkingI 19841 |Lada|1987| ) and to be the precursors of the de- 



* Based on observations made with ESO Telescopes at the La Silla 
or Paranal Observatories under programme ID 380.C-0326. 
** e-mail: cgraefe@astrophysik.uni-kiel.de 



let al.|2009| ). 



Investigating the evolution and thus dissipation of circum- 
stellar disks is crucial for our understanding of planetary system 
formation. Various efforts to spatially resolve the inner struc- 
ture of transitional disks have been undertaken (e.g. with the 
Submillimeter Array by [Brown et al.||2009] ). However, owing 
to their limited sensitivity, submillimeter and millimeter obser- 
vations only provide information about the large reservoirs of 
cold dust in the outer regions so far. Thus, high-resolution mid- 
infrared imaging of transitional disks is essential to obtain infor- 
mation about the warm dust in the innermost regions of the disk, 
to detect inner-disk clearings as signatures of the disk dissipation 
mechanisms and to investigate the structure close to the star and 
therefore the regions where planets potentially form (e.g., Wolf] 
let al.|2007D . 
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Table 1. Observing log 



Target 


RA (J2000) 


Dec (J2000) 


Observation Date 


Observation Time 


Exposure Time 




[hli mm ss] 


[dd mm ss] 




[hh:mm] 


[s] 


DHTau 


04 29 41.56 


+26 32 58.3 


2007 Oct 24 


07:27 


883.2 


DMTau 


04 33 48.72 


+18 10 10.0 


2007 Oct 20 


06:26/07:30 


2x1766.4 


GMAur 


04 55 10.98 


+30 21 59.5 


2007 Oct 20 


08:30 


353.3 



Table 2. Properties of the observed standard-stars 



Target 


RA (J2000) [hh mm ss] 


Dec (J2000) [dd mm ss] 


Spec. Type 


msiv-Fiiter [mag] 


HD31421 


04 56 22.27 


+13 30 52.1 


K2III 


11.38 


HD 27482 


04 21 15.26 


+27 21 00.9 


K5III 


6.12 



Notes. From flux catalog of mid-infrared standard-stars for VISIR imager filters based on |Cohen et al.| ( [T999) . 



In the mid-infrared, various efforts to spatially resolve cir- 
cumstellar disks have been undertaken. For example, Chen et al. 
([2008) and Churcher et al.| (| 2011| presented resolved mid- 



infrared images of the disk of HD 191089 and HD 181327, re- 
spectively. The circumstellar disk of the well-studied Herbig 
Ae/Be star AB Aur was investigated by |Chen & Jura| ( |2003| ) and 
[Pantin et al.| p005 ) in the mid-infrared. 

In contrast to previous mid-infrared imaging of circumstellar 
disks, our observations and data analysis presented here reach 
the limit of angular resolution of today's single-dish, ground- 
based telescopes. In the case of AB Aur, which emits a N band 
flux of >20 Jy, structures with a radius of ~280 AU were detected 
( [Chen & Jura|2003[ [Pantin et al.|2005| ). We investigate sources 
that are -100 times fainter and structures that have more than 
ten times smaller radii. Although we are using one of the best ob- 
serving facilities available, the properties of the observed objects 
necessitate a new data analysis procedure, which is described in 
Sect. 4. An even higher resolution in this wavelength range can 
only be achieved by means of interferometric observations, e.g., 
with MIDIA^LTI (e.g. Schegerer et al.'2009 ). However, such ob- 
servations are still limited to the mid-infrared brightest circum- 
stellar disks. 



2. Target objects 

Our target objects are three transitional disks that were selected 
based on their 10 ^m flux from a list of most well-studied poten- 
tial transitional di sks ( [Forrest et al.|2004t[Hartmann et al.|2006t 
Furlan et al. 2006 1. These three objects are DH Tau, DM Tau, and 
GIVI Aur, which are all located in the Taurus/Auriga star- forming 
region at a distance of ~140pc (Kenyon et al. 1994). Their SEDs 
(see Fig. 1 in [Furlan et al.|2006j show significant infrared excess 
emission, which is at least ten times higher than the stellar emis- 
sion in the mid-infrared. As it is typical of transitional disks, 
their SEDs also have a flux deficiency in the near-infrared with 
increasing slopes in the mid-infrared wavelength range, indicat- 
ing an inner radius of the dust disk much larger than the dust 
sublimation radius. 

In the following, our target stars are briefly characterized: 



1 . DH Tau is a Class II source of spectral type M 1 ( Watson et al. 



2009 1. Itoh et al. (2005 1 detected a young brown dwarf com- 
panion to the transitional disk around DH Tau, separated by 
2.3", corresponding to a spatial distance of 330 AU if the 
orbit is seen face-on. 



DM Tau is a low-mass, pre-main sequence star that does not 
show any significant excess emission below 8 \im, but is sur- 
rounded by a still accreting transitional disk with a mass ac- 
cretion rate slightly lower than t he average of disks in Taurus 



average 



= 10"^-10"^MQ/yr) (iHartmann et al. 



1998 



White 



f& Ghez 2(301] [SaTykeTal. 2009). Its inner-dis k hole radius 



is expected to be 3 AU with little dust inside (Calvet et al 
.2005, Pott etal |2010| ). 

GM Aur is a 3 Myr old K7 T Tauri star that is surrounded 
by a transition al disk and is also known to be still in its ac- 
creting phase ( [Calvet et aL[[2005[ [Dutrey et aL][2008 |). Its 
SED indicates the presence of an inner-disk hole of radius 
20 - 24 AU that is partially filled with optically thin dust 



(Calvet et al. 



2005 Hughes et al. 2009 Pott et al. 2010 1 



Dutrey et aLf( 2008) and [Hughes et al.[( 2009|l suggest that 



the inner-disk hole plausibly results from the dynamical in- 
fluence of a planet on the disk material. 

3. Observations and data reduction 

The mid-infrared observations of our three targets were carried 
out on 2007 October 20 and 24, as described in Table [T[ us- 
ing VISIR ( [Lagage et al. 1200 4) in the SIV filter (A, = 10.49 \xm, 
A/l = 0.16 [xm). VISIR is the ESO/VLT mid-infrared imager and 
spectrograph, composed of an imager and a long-slit spectrom- 
eter covering several filters in N and Q bands. To achieve the 
highest possible angular resolution, the small field objective was 
used (0.075'7pixel). A standard chop/nod-procedure to remove 
sky and telescope emission with a chop of 8" in the north-south 
direction and a perpendicular 8" nod were applied for the ob- 
servations. Calibration observations of standard- stars within a 
few degrees of the science object were taken immediately be- 
fore the science observations. The standard-stars were chosen 
from the mid-infrared spectrophotometric standard-star catalog 
of the VLT based on [Cohen et aL] ( [T999[ ) (see Table [2]). These 

Table 3. Observation properties of the observed standard- stars 



Target 


Obs. Date 


Obs. Time 


Exposure 


NBand 




[2007] 


[hh:mm] 


Time [s] 


Seeing ["] 


HD31421 


Oct 20 


06:16 


61.44 


0.29 


HD 31421 


Oct 20 


07:22 


61.44 


0.26 


HD 27482 


Oct 20 


08:20 


61.44 


0.29 


HD 27482 


Oct 24 


07:16 


61.44 


0.29 
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standard-stars were used to characterize the point spread func- 
tion (PSF) for comparison with the science sources to detect any 
signs of a circumstellar disk. All observations were executed in 
service mode. Raw data were reduced using the VISIR data re- 
duction pipeline provided by ESO. A full-resolution background 
map was created f or every observed ob ject using the Source 
Extractor software (Bertin & Amouts"1996). This map was sub- 
tracted from the image to remove the background. This approach 
provides a more reliable estimate of the background than the as- 
sumption of a constant background. 

4. Data analysis 

4. 1. Radial flux profile 

The key components of each of our science targets are the cen- 
tral star and the surrounding disk. In contrast to the brightness 
profile of the corresponding reference (standard-) star, the pro- 
files of the science targets are extended by the circumstellar disk. 
The standard-stars, which are considered to be point sources, are 
used as PSF references (see Tables [2] and [3]). To obtain the flux 
profile of a science object and a corresponding standard- star, we 
carry out a radial averaging of each object. For this purpose, 
the images are resized by a factor of ten using bilinear inter- 
polation. The center of each object is determined by fitting a 
two-dimensional Gaussian profile. The radial averaging is then 
performed with respect to the obtained object center and yields 
as a result the radial flux profile of the object. Assuming radial 
symmetry, this averaging method yields a much better flux pro- 
file than only one slice through the object and also considerably 
improves the signal-to-noise ratio (hereafter S/N) by a factor of 
ViV, where is the number of pixels in every annulus. Radially 
symmetric disk images can of course only be expected if the 
disks are seen in face-on orientation. The efl'ect of a disk incli- 
nation other than face-on on our results is discussed in Sect. 4.3. 

Owing to the small chopping and nodding throw in the ob- 
servations, the observed source is always present in the field of 
view of the instrument such that the source is detected four times 
in the final reduced image (as two positive and two negative im- 
ages). To increase the S/N by a factor of two, the radial averaging 
is accomplished for all four observations of the object in the final 
image and averaged afterwards. This profile is then normalized 
to a peak flux of 1. Subtracting the standard-star profile from 
the corresponding science object profile yields the flux residuals 
produced by the extended star-disk system. 

4.2. Detection threshold 

At the next step, one has to investigate whether the derived flux 
residuals are significant. Therefore, we define a detection thresh- 
old such that any residual flux above this threshold is caused pri- 
marily by an extension of the star-disk system, i.e., a spatially 
resolved circumstellar disk surrounding the star. 

Our detection threshold takes into account both the standard 
deviation in the background flux (cTbg) of the science image and 
the seeing variations during the observations of the standard- 
star and science object. The background noise is calculated from 
the full-resolution background noise map created using Source 
Extractor. We impose the constraint that the residual flux must be 
higher than Scr^g. Owing to the radial averaging, the background 
noise of the flux profile depends on the radius and declines with 
increasing radius. 

Standard-stars are always observed before and/or after the 
science object to have nearly the same observing conditions, thus 



for point sources nearly the same fuU-width-at-half-maximum 
(FWHM). However, the FWHM depends on the seeing, which 
varies during observations. If not properly taken into account, 
variations in seeing can cause a difl"erence signal similar to the 
case of an extended star-disk system, thus mimic a spatially 
resolved disk. For our data analysis, it is therefore crucial to 
determine the seeing variations during the observations of the 
standard-star and science object. This is difficult because infor- 
mation about the seeing during observations of the science ob- 
ject at the observed wavelength is unavailable. Only the FWHM 
of the standard- star is known from which the seeing at both the 
observation time and the observed wavelength must then be de- 
termined (see Table[3]). To estimate the seeing variations, we an- 
alyze observations of our standard-starHD 31421, which was re- 
peatedly observed over several hours before and after the obser- 
vations of a scientific target. Therefore, we use the observation 
with the highest and the smallest FWHM of this standard- star, 
i.e., at the worst and the best seeing conditions and determine 
the radial flux profile of both as described above. Assuming that 
the seeing has a normal distribution during the observations, the 
seeing variations can now be calculated as the half difl^erence 
of these two profiles. This method yields a radial profile reflect- 
ing a conservative estimate of the seeing variations during the 
observations. The derived seeing variations are adopted for all 
observed science objects. 

To increase the statistical significance of the derived flux 
residuals, we also insist that they be three times larger than the 
derived seeing variations. Our detection threshold is now defined 
as the sum of 3crbg and three times the seeing variations. The 
contribution of the two components to this threshold is exemplar- 
ily represented for GM Aur in Fig.[T] where it is clearly visible 
that the seeing variations are a significant part of the threshold, 
hence must be reliably determined. 

4.3. Model 

We present a simple model to determine the inner-disk hole ra- 
dius of our objects, provided that they are spatially resolved in 
our observations. However, our approach also allows us to esti- 
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Fig. 1. The contribution of the two factors to the detection thresh- 
old. 
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mate an upper limit to the inner-disk hole radius of our transi- 
tional disks that could not be spatially resolved. 

Our model consists of a star surrounded by a disk observed at 
a certain inclination. To allow comparison with the observations, 
the image of this model was convolved with the measured PSF. 
This PSF is the arithmetic mean of the four observations of the 
observed standard-star in the final image, which are caused by 
the chop/nod technique. For every model of the science objects, 
the corresponding standard-star is used. 

The inner disk rim, which dominates the mid-infrared sur- 
face brightness distribution of the disk, is represented by an in- 
clined circle with a radius equal to the gap/hole radius of the 
transitional disk whose inclination equals the known inclination 
of the observed disk. The star is located in the center with a 
star/disk flux ratio obtained from the observed SED ( [Furlan et al.| 
2056} [Hughe s et al. 2009) (see Table |4]). As described in Sect. 
4.1, the flux profiles of the PSF-convolved object and the mea- 
sured PSF are obtained using a radial averaging method. They 
are normalized to a peak flux of 1 . The radial profile of the mea- 
sured PSF is then subtracted from the convolution profile to ob- 
tain the flux residuals arising from the extended star-disk sys- 
tem. These flux residuals and the determined detection threshold 
yield the model to the corresponding observation. 



Table 4. Star/disk flux ratio at 10 [xm of the observed objects. 
For details, see Sect. 4.3. 



Table 5. Constraints on the inner-disk hole radii and photometry 
of the observed transitional disks 



Target 


Star/disk flux ratio 


DHTau 


0.111±0.018 


DMTau 


0.125±0.015 


GMAur 


0.081±0.013 



4.4. Photometry 

The N band fluxes of the three observed objects are determined 
using the VISIR photometry package by [Honig et al. (2010 1. 
For photometry, the science targets can be considered as point 
sources such that each science and calibrator beam is fitted by a 
two-dimensional Gaussian. The fluxes are then obtained by cal- 
culating conversion factors from the integrated intensity of the 
Gaussian fits of the calibrators and using these factors to eval- 
uate the respective integrated Gaussian intensity of the science 
targets. 



5. Results and discussion 

The resulting constraints on the inner-disk hole radii and the 
photometry of the observed transitional disks are compiled in 
Table [5] In the following, we discuss the analysis of the individ- 
ual objects in detail. 

For GM Aur, we are able to detect the extension of the star- 
disk system caused by its circumstellar disk (see Fig. [2j left). 
Estimates of the disk inclination of GM Aur range from 49° to 
56° (|Simon e t al. 2000; Calvet et al. 2005; Du trey et al.|[2008t 
[Hughes et al.t2009,,Isella et al. 2 009). Comparing the observa- 
tion with our model, we can constrain the inner-disk hole radius 
of GMAur's transitional disk to 20.5+^° AU (see Fig. [2j right), 
where the smaller value represents the radius for the smaller in- 
clination and the larger value the radius for the larger inclina- 
tion. This size is comparable to the observational seeing in the 



Target 


Inner-disk hole radius [AU] 


Flux [mJy] 


DHTau 
DMTau 
GMAur 


<15 5+°'5 
20 5"^' ° 


178 ±31 
56 ±6 
229 ± 14 



Notes. The uncertainties in the inner-disk hole radii are caused mainly 
by the uncertainty in the disk inclination (see Sect. 5). 



N band as a disk hole of radius 20.5 AU at 140 pc amounts to 
~ 0.29" across. The inner-disk hole radius derived from our ob- 
servations is in very good agreement with the values of ~20 AU 
and ~ 24 AU found by [Hughes et al] ( [2009 ) and C alvet et al. 
( [2005 1 , respectively. 



In Figs. [3] and [4| we present the results and models for 
DHTau and DMTau. For both objects, the difference in the ob- 
ject and standard-star flux profiles is below the detection thresh- 
old, i.e., no significant extension of the star-disk system is de- 
tected. Using our model, we can investigate at which inner-disk 
hole radius we would have detected a significant extension of the 
star-disk system, hence measure an upper limit to the inner-disk 
hole radius. The disk inclination of DMTau ranges from 32° to 
40° (Simon et al. 2000; Calv et et al..2005,|Pott et al.|2010| ). We 
therefore infer that DMTau's inner-disk hole radius is smaller 
than 15.5;^QgAU, which thus agrees \ 
who expected an inner-disk hole radius of 3 



Calvet et al. 
AU.ForUH 



(20051 



lau, no 

estimate of the disk inclination has yet been made. Hence, using 
our model we determine an upper limit to the disk hole radius 
of 15.5;|;2 o^U, where the lower limit corresponds to when the 
disk is seen face-on and the upper limit to an edge-on orienta- 
tion. However, a close to edge-on orientation of this disk can 
be excluded because the mid-infrared bright inner-disk region 
would not then be visible through the optically thick disk. No 
complementary constraints on the inner-disk hole radius are yet 
available for DH Tau. 

When considering deviations in both the radial flux profiles 
and the detection threshold, besides the detection threshold per 
se, the error bars are far too small to be visible in the plots. The 
error in the radial profiles is the background noise, which is im- 
plemented in the detection threshold. The error in the second part 
of this threshold, i.e. in the seeing variations, is half times the 
sum of the quantity crbg of the standard-star for the worst seeing 
conditions plus the one for the best. Owing to their high S/N and 
the radial averaging, the background noise of the standard- stars 
is negligibly small. Therefore, this also applies to the error in the 
seeing variations. The contribution of 3crbg and three times the 
seeing variations to the detection threshold can be seen in Fig.[T] 

The performed N band photometry of our three observed 
transitional disk objects yields fluxes of 178 + 31mly for 
DH Tau, 56 + 6 miy for DM Tau, and 229 + 14 miy for GM Aur. 
They are the first imaging photometric measurements of these 
objects at this wavelength. As a complement to already known 
fluxes at other wavelengths, they will be useful to constrain the 
SEDs of these three objects and model their circumstellar disks. 

Spatially resolved 2 \.im Keck Interferometer observations of 
DM Tau and GM Aur presented by [PottetaL] ( [2010| ) indicate that 
there is a significant K band flux contribution from extended hot 
circumstellar matter very close to the star. This is consistent with 



previous SED modeling of DMTau and GMAur ( Calvet et al 
2005 ; Hu ghes et al.|2009] ) that also inferred that the amount of 
dust in the inner few AU is too small to contribute significantly 
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Fig. 2. Left: Observed flux residuals of GM Aur including the determined detection threshold. Right: Modeled flux residuals for 
several inner-disk hole radii including the detection threshold. 
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Fig. 3. Left: Observed flux residuals of DH Tau including the determined detection threshold. Right: Modeled flux residuals for 
several inner-disk hole radii including the detection threshold. 



to the mid-infrared flux. Nevertheless, our results help to reduce 
the ambiguities of the SED modeling and allow us to conclude 
that this hot dust has to be located inside our measured inner- 
disk hole radii. Higher-resolution observations, which will po- 
tentially be feasible with the Atacama Large Millimeter Array, 
will allow us to spatially resolve substructures within ~ 10 AU 
and thus investigate whether the inner thin disk traced in K band 
is connected to the "outer" disk or a gap exists. 
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Fig. 4. Left: Observed flux residuals of DM Tau including the determined detection threshold. Right: Modeled flux residuals for 
several inner-disk hole radii including the detection threshold. 
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